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Abstract 



Weak localization leads to the same correction to both the conductivity 
and the McMillan's electron-phonon coupling constant A (and Xtr)- Conse- 
quently the temperature dependence of the thermal electrical resistivity is 
decreasing as the conductivity is decreasing due to weak localization, which 
results in the decrease of the temperature coefficient of resistivity (TCR) with 
increasing the residual resistivity. When A and Xt r are approaching zero, only 
residual resistivity part remains and gives rise to the negative TCR. Accord- 
ingly, the Mooij rule is a manifestation of weak localization correction to the 
conductivity and the electron-phonon interaction. This study may provide a 
new means of probing the phonon- mechanism in exotic superconductors. 



PACS numbers: 72.10.Di, 72.15.Rn, 72.15. Cz, 72.60.+g 
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I. INTRODUCTION 



Although weak localization has greatly deepened our understanding of the normal state 
of disordered metals, 1 ' 2,3 its effect on superconductivity and the electron-phonon interaction 
has not been understood well. 2 Recently, it has been shown that weak localization leads 
to the same correction to the conductivity and the phonon- mediated interaction. 4 ' 5 In fact, 
there are overwhelming numbers of experiments which support this idea. 5 For instance, 
tunneling, 6-8 specific heat, 9 x-ray photoemission spectroscopy (XPS), 10 correlation of T c 
and the residual resistivity, 11-13 universal correlation of T c and the resistance ratio, 14-16 
and loss of the thermal electrical resistivity 17 with decreasing T c clearly show a decrease 
of the electron-phonon interaction accompanying the decrease of T c with disorder. It is 
then anticipated that the electron-phonon interaction in the normal state of metals will 
also be influenced strongly by weak localization. We expect that phonon-limited electrical 
resistance, attenuation of a sound wave, thermal resistance, and a shift in phonon frequencies 
may change due to weak localization. 18 

Indeed, the Mooij rule 19 in strongly disordered metallic systems seems to be a manifesta- 
tion of the effect of weak localization on the electron-phonon interaction and the conductivity. 
In early seventies, Mooij found a correlation between the residual resistivity and the tem- 
perature coefficient of resistivity (TCR). In particular, TCR is decreasing with increasing 
the residual resistivity. Then it becomes negative for resistivities above 150/ifkm. We stress 
that this behavior is consistent with the above superconducting properties: correlation of 
T c and the residual resistivity, 11-13 universal correlation of T c and the resistance ratio, 14-16 
and loss of the thermal electrical resistivity 17 with decreasing T c . 

There are already several theoretical works on this problem. Jonson and Girvin 20 per- 
formed numerical calculations for an Anderson model on a Cayley tree and found that the 
adiabatic phonon approximation breaks down in the high- resistivity regime producing the 
negative TCR. Imry 21 pointed out the importance of incipient Anderson localization (weak 
localization) in the resistivities of highly disordered metals. He argued that when the inelas- 
tic mean free path, £ p h, is smaller than the coherence length, £, the conductivity increases 
with temperature like and thereby leads to the negative TCR. On the other hand, Kaveh 
and Mott 22 generalized the Mooij rule. Their results are as follows: The temperature depen- 
dence of the conductivity of a disordered metal as a function of temperature changes slope 
due to weak localization effects, and if interaction effects are included, the conductivity 
changes its slope three times. Gotze, Belitz, and Schirmacher 23 ' 24 introduced a theory with 
phonon-induced tunneling. There is also the extended Ziman theory, 25 and Jayannavar and 
Kumar 26 suggested that the Mooij rule can arise from strong electron-phonon interaction 
taking into account qualitatively different roles of the diagonal and off-diagonal modulations. 

In this paper, we propose an explanation of the Mooij rule based on the effect of weak 
localization on the electron-phonon interaction. If we assume the decrease of the electron- 
phonon interaction due to weak localization, 4,5 we can understand the decrease of TCR with 
increasing the residual resistivity. The negative TCR is therefore due to weak localization 
correction to the Boltzmann conductivity, since when TCR is approaching zero there is 
no temperature-dependent resistivity left. (This latter point is similar to Kaveh and Mott's 
interpretation. 22 ) In Sec. II, we briefly describe the Mooij rule. In Sec. Ill, weak localization 
correction to the McMillan's electron-phonon coupling constant A and \ tr is calculated. A 
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possible explanation of the Mooij rule is given in Sec. IV, and its implication is briefly 
discussed in Sec. V. In particular, this study may provide a means to probe the phonon- 
mechanism in exotic superconductors. 

II. THE MOOIJ RULE 

Mooij 19 was the first to point out that the size and sign of the temperature coefficient 
of resistivity (TCR) in many disordered systems correlate with its residual resistivity po as 
follows: 

dp/dT > if po < Pm 

dp/dT<0 if p >pu- (1) 

Thus, TCR changes sign when p reaches the Mooij resistivity pu — l50pQcm. An approx- 
imate equation for p(T) is given by 2 

p(T)=p + (p M -p )AT, (2) 

where A is a constant which depends on the material. 

Figure 1 shows the temperature coefficient of resistance a versus resistivity for transition- 
metal alloys obtained by Mooij. It is clear that a (and TCR) is correlated with the residual 
resistivity. Note that above 150pQcm most a's are negative while no negative a is found for 
resistivities below 100pf2cm. Figure 2 shows the resistivity as a function of temperature for 
pure Ti and TiAl alloys containing 3, 6, 11, and 33% Al. TCR is decreasing as the residual 
resistivity is increasing. For TiAl alloy with 33% Al shows the negative TCR. We note that 
the positive TCR is basically high temperature phenomenon, presumably related to the 
phonon-limited resistivity, whereas the negative TCR is rather low temperature behavior, 
probably connected with the residual resistivity part. Since this behavior is generally found 
in strongly disordered metals and alloys, amorphous metals, and metallic glasses, 2 it is called 
the Mooij rule. However, the physical origin of this rule has remained unexplained until now. 

III. WEAK LOCALIZATION CORRECTION TO THE ELECTRON-PHONON 

INTERACTION 

Since the electron-phonon interaction in metals gives rise to both the (high temperature) 
resistivity and superconductivity, these properties are closely related, which was noticed by 
many workers. 27-31 Gladstone, Jensen, and Schrieffer 27 pointed out that A and the high 
temperature electrical resistivity are closely related each other. Hopfield 28 ' 29 noted that 
the electronic relaxation time due to electron-phonon interaction, as measured in optical 
experiments above the Debye temperature, should be approximately equal to 2n\kBT/h. 
He applied this idea to Nb, Mo, Al and Sn and found a good agreement with experiment. 
Grimvall 30 estimated A for noble metals from Ziman's high temperature resistivity formula. 
Maksimov and Motulevich 31 followed the idea of Hopfield and estimated A from optical 
measurements for Pb, Sn, In, Al, Zn, Nb, V, NbsSn, and "V^Ga, which are in good agreement 
with the McMillan's coupling constant A from superconductivity data. 

In this Section, we show that weak localization leads to the same correction to the 
conductivity, the McMillan's electron-phonon coupling constant A, and X tr . 
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A. High Temperature resistivity 



At high temperatures, the phonon limited electrical resistivity is 



32-35 



Anmk B T f a 2 tr F(u) 



Pph (T) = ^- / ^-^dco, 

_ 2ixmk B T 

— 57 — Mr, [o) 

where a tr includes an average of a geometrical factor 1 — cos9^, and F(u) is the phonon 
density of states. On the other hand, in the strong-coupling theory of superconductivity, 36,37 
the McMillan's electron-phonon coupling constant is defined by 37 

A = 2 / — v ; v ' duo. (4) 
J ui 

Assuming o? tr = a 2 , 32 ' 38-40 we obtain 

2nmk B T 

PAT) = -^Ij-^ (5) 
= 27rmk 2 f X . (6) 

Consequently the McMillan's coupling constant A also determines the size and sign of TCR. 

The existence of this relationship was well confirmed theoretically and experimentally. 
Table I shows the comparison of X tr and A by Economou 38 for various materials. He obtained 
X tr from Eq. (5) and compared with A, as obtained from T c measurements, and/or tunneling 
experiments, and/or first principle calculations. 39 The overall agreement between \ tr and 
A is impressive. Grimvall estimated A for noble metals 30 and noble metal alloys 41 from 
Eq. (6). Maksimov 40 also noted the direct relation between A and the high temperature 
resistivity. Hayman and Carbotte 42 pointed out that information on the volume dependence 
of an electron-phonon coupling strength can be obtained from high temperature resistivity. 
Chakraborty, Pickett, and Allen 43 used Eq. (5) to obtain the empirical values of A tr for Nb, 
Mo, Ta, and W. They found that \ tr from resistance and the McMillan's coupling constant 
A from superconductivity are very similar in magnitude for these materials. We can also 
mention experimental confirmations by Rapp and Crawfoord 44 for Nb-V alloys, Rapp and 
Fogelholm 45 for Al-Mg alloys, Fliikiger and Ishikawa 46 for Zr-Nb-Mo alloys, Fogelholm and 
Rapp 47 for In-Sn alloys, Lutz et al. 48 for Nb 3 Ge films, Man'kovskii et al. 49 for thin Sn films, 
Rapp, Mota and Hoyt 50 for Au-Ga alloys, and Sundqvist and Rapp 51 for aluminum under 
pressure. Figure 3 shows the McMillan's coupling constant A versus dp/dT oc Xt r for Au-Ga, 
Au-Al, and Ag-Ga alloys, 52 which exemplifies the correlation implied by Eq. (6). 



B. Weak localization correction to the McMillan's coupling constant A and X tr 

Now we need to calculate the McMillan's electron-phonon coupling constant A for highly 
disordered systems. We follow McMillan's approach to the strong-coupling theory. 37 ' 5 (For 
simplicity we consider an Einstein model with frequency uj d ). He showed that A can be 
written as 37 
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x=2f^Mm duJ (7) 

J Id 

< I 2 > 

where M is the ionic mass and iVo is the electron density of states at the Fermi level. 
< I 2 > is the average over the Fermi surface of the square of the electronic matrix element 
and < u 2 >= oj 2 d . In the presence of impurities, weak localization leads to a correction to 
a 2 (uo) or < I 2 >, (disregarding the changes of F(uo) and N ). 

There are two ways to obtain the McMillan's coupling constant A in the presence of 
impurities. One method is to calculate A directly from Eq. (8), using the electronic matrix 
element for disordered systems and the other is to carry out the canonical transformation of 
Frohlich in the scattered state basis. 5,53 We have found that both methods lead to the same 
A. 

In this paper, we use the latter method in a simple manner by observing that the Frohlich 
interaction can be derived from the phonon Green's function. 54 We note that the equivalent 
electron-electron potential in the electron-phonon problem is given by the phonon Green's 
function D(x - a;'): 54 " 56 

Viz-aO^j^Dix-aO, (9) 

where x = (r, t) and Iq is the electronic matrix element for the plane wave states. The 
Frohlich interaction at finite temperatures is then given by 54 

V nn ,(u,u') = j^Lr J J dvdv'r n ^WnV)D{v - r> - w> fi (r>„(r) 

7~2 2 
/ I /. /"_M2|„/. /„\|2 t U D 



Mu 2 D J ^' (r)l '^ (r)|rfr 4 + (---0 2 

UJjj + (id — Id' 



Vnn' o , " ,, A2 , (10) 



where 



D{v-v\d-u')=Y. 



Up c iq-(r-r>) 



(d 



Ld') 2 +id 2 D 



id 2 



(ou - d') 2 + u 2 D 



Here u means the Matsubara frequency and ip n and vfj n denote the scattered state and its 
time-reversed partner, respectively. Therefore, we get the strong-coupling gap equation 5 

1 2^1 (, , , ,A2 i , ,2 v nri , 2 , p 2 /, ,/\ ' ^ 1Z V 
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where 

EnW) = y/el + AKu'), (13) 
and the McMillan's electron-phonon coupling constant A 

A = N < V nn ,(0,0) >= N j^- < J |Vn(r)|Vn'(r)| 2 rfr > . (14) 

Here e n means the eigenenergy of the scattered state tp n . As expected, Eq. (14) is the 
same as Eq. (8). It is remarkable that the McMillan's electron-phonon coupling constant is 
determined by the density correlatin function. 

Note also that in the presence of impurities, the density correlation function has a free- 
particle form for t < r (scattering time) and a diffusive form for t > r. 57 As a result, for 
t > t (or r > £), one finds 58-62 

R(t>r)= I |V>„(r)|> n ,(r)| 2 dr 

Jt>T 

= Y / \<*Pn\e i niJn>>\ 2 

q 

= 2*WD? (15) 

(i-t)- ( 16 ) 



2(k F £) 2K V 

Here I is the mean free path and L is the inelastic diffusion length. Whereas the contribution 
from the free-particle-like density correlation for t < r is 5 ' 58 

R(t<r) = ( |^ n (r)|V n ,(r)| 2 dr 

Jt<T 

(17) 

Since the phonon-mediated interaction is retarded for t ret ~ 1/wd, only the free-particle-like 
density correlation contributes to A. This is also true of X tr , simply because the conductivity 
is determined by the behavior of the wavefunction ip for t < r (or r < £). m 

Consequently, we obtain weak localization correction to the McMillan's coupling constant 
A and \ tr - 

^^-(wOi- < 18 > 

and 



Muj 2 d 1 (k F £) 2 y V 



N I£ r . 3 



We have used the fact that L is effectively infinite at T = 0. Note that the weak localization 
correction term is the same as that of the conductivity. 
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IV. EXPLANATION OF THE MOOIJ RULE 



As noted in the Section II, the positive TCR is high temperature phenomenon whereas 
the negative TCR is low temperature phenomenon. Thus, the decrease of the positive TCR 
is mainly due to the decrease of the phonon-limited resistivity, which is a manifestation of 
weak localization correction to the electron-phonon interaction. On the other hand, the 
negative TCR originates from the residual resistivity, which is also a manifestation of weak 
localization correction to the conductivity. Accordingly, weak localization seems to be the 
physical origin of the Mooij rule in disordered metals. One should note that this observation 
agrees with the superconducting behavior of disordered systems, such as correlation of T c 
and the residual resistivity, 11-13 universal correlation of T c and the resistance ratio, 14-16 and 
loss of the thermal resistivity 17 with decreasing T c . 



A. Decrease of TCR at high temperatures 

Upon substituting Eq. (19) into Eq. (3), one finds the phonon-limited high temperature 
resistivity 

2ixmk B T 

^ 2nmk B TN I 2 3 

ne 2 h Mcu 2 D [ (k F £) 21 ' { ' 

Note that as the disorder parameter l/k F £ is increasing, both the magnitude of the phonon- 
limited resistivity and the TCR decrease. This behavior is due to the reduction of the 
McMillan's electron-phonon coupling constant when electrons are weakly localized. It is 
remarkable that the slope of the high temperature resistivity varies as ~ l/(k F £) 2 , in accord 
with the behavior of the residual resistivity. 

The phonon-limited resistivity p p h versus temperature T is shown in Fig. 4 (a) for six 
values of k F £. We used k F = 0.8l -1 , n = k 3 F /3Tc 2 , and N I^/(Muj 2 d ) = 0.5. It is clear that 
TCR is decreasing significantly as the electrons are weakly localized. 



B. Negative TCR at low temperatures 

At low temperatures the conductivity and the residual resistivity are given by 2 ' 3 

^i 1 -^ 1 -!"- < 2i > 

and 

Po = ~T7i 3 n FTT' ( 22 ) 

where a B = ne 2 T/m. When l/k F £ becomes comparable to ~ 1, the magnitude and slope of 
pphiT) are negligible. In that case, only the residual resistivity will play an important role. 
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Therefore, the observed negative TCR may be understood from the residual part. With 
decreasing T, since the inelastic diffusion length L increases, the residual resistivity will also 
increase, leading to the negative TCR. We stress that both the phonon-limited resistivity 
and the residual resistivity have the same quadratic dependence on the disorder parameter 
l/k F £. 

Figure 4 (b) shows the temperature dependence of the residual resistivity po for kpi = 
2.2, 2.4, 2.8, 3.4, 5, and 15. Since it is difficult to evaluate kpi up to a factor of 2, 64 we 
assumed that p = lOOpQcm corresponds to kpi = 3.2. We used the same kp as in Fig. 4 
(a) and L = \ZDt,i = \/i x 350/T(A). Here D is the diffusion constant and T; denotes the 
inelastic scattering time. When kpi is comparable to 1, the negative TCR emerges. Notice 
the scale difference between Figures, 4 (a) and 4 (b). 



C. Comparison with experiment 

In Sections IV. A and B, we have explained the physical origin of the Mooij rule. In this 
section, we compare our theoretical resistivity curve and the experimental data (Figure 2) 
for extended temperature range. Let us remind the approximate formula for p(T) suggested 
by Lee and Ramakrishnan, i.e., 2 

p(T) =p + (p M - Po)AT. (23) 

This form of equation can be obtained if we add the residual resistivity Eq. (22) and the 
phonon-limited resistivity Eq. (20), that is, 



p(T)= Po + Pph (T) 



^[l-O^(l-f)] w*h M< (M) 2J - 

It should be noticed that the addition of both resistivities does not mean the Matthiessen's 
rule. Here we included the interference effect between the electron-phonon and electron- 
impurity interactions: 

p(T) = Po + p ph (T, c = 0) + A$f , (25) 

where c denotes an impurity concentration. Whereas Altshuler 65 and Reizer and Sergeev 66 
investigated corrections to the impurity resistivity due to the interference, we have considered 
its correction to the phonon-limited resistivity. Since the interference correction to the 
impurity resistivity is ~ 1% of the residual resistivity, 66 ' 67 we neglect its effect for simplicity. 
In general, the phonon-limited resistivity at any temperature T is given by 32-35 

PAT) = ^ / mu)aUu)F^) 
Fp y ' ne 2 J (el 3 *" - - e-P*") 

where f3 = 1/ksT. For an Einstein phonon model with 68 

at n 

a£HFH = ^%-^), (27) 
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it is rewritten as 69 

27rm jV I 2 ({3hu D )m D 
Pph[ ' ne 2 Mo; 2 , (e^ - 1)(1 - e~^) ' 1 J 

It is necessary to emphasize that this result is exact for the phonon-limited resistivity in an 
Einstein model. 67 Including the weak localization correction to a 2 (uo) = a 2 r (u), 

^™ = ^ 1 -<W Wu - UD) - (29) 

one finds 

_ 2?rm jV / 2 3 (phu D )u D 

Pph{ ' ne 2 Muj 2 d [ (k F £) 21 (e^ - 1)(1 - e~^) ' { } 

Finally, we obtain the total resistivity at any temperature T: 



p(T)=p + p ph (T) 



1 2nmNoll 3 (phcu D )u D 



*b[1- (j^(l-f)] ne 2 M^ L (M) 2J (e^ - 1)(1 - e -^) ' 

(If we consider the Debye and realistic phonon models, there are minor changes. However, 
the overall behavior is the same. More details will be published elsewhere.) 

Figure 5 shows the resistivity as a function of temperature for k F £ = 2.3, 2.5, 2.8, 3.4, 5, 
and 15. The solid lines represent the resistivity from an accurate expression Eq. (31), while 
the dashed lines are obtained from Eq. (24). We used the same parameters as those in Fig. 
4 and huo = 250fT. It is noteworthy that both equations give rise to almost the same curve 
as the system is more disordered. For low temperatures is determined by electron-electron 
scattering while for high temperatures it is determined by the electron-phonon scattering. 
Since we are interested in rather high temperatures, we assumed r, ~ T~ l corresponding to 
the electron-phonon scattering. 2,3 Considering the crudeness of our calculation, the overall 
behavior is in good agreement with experiment, Fig. 2. 



V. DISCUSSION 

At low temperatures the interference of the Coulomb interaction and the impurity scat- 
tering leads to the interaction correction to the conductivity 60,2 This effect is described 
by™ 

a =^ i -v^ {i --L ) -(£ef {i -£ )l (32) 

where L T = (hD jk B T) x l 2 and C ~ 1. The second correction term is the interaction term. 
The constant C, however, changes sign depending on the exchange and Hartree terms and 
since it is difficult to determine C, 2 ' 3,70 we did not include this term. But it may be important 
at much lower temperatures. 
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It is clear that weak localization effect on the electron-phonon interaction needs more 
theoretical and experimental studies. In particular, weak localization effect on the attenua- 
tion of a sound wave, shear modulus, thermal resistance, and a shift in phonon frequencies 
will be very interesting. Since superconductivity is also caused by the electron-phonon in- 
teraction, comparative study of the normal and superconducting properties of the metallic 
samples will be beneficial. There is already compelling evidence that this is the case. For 
instance, Testardi and his coworkers 14-17 found the universal correlation of T c and the re- 
sistance ratio. They also found that decreasing T c is accompanied by the decrease of the 
thermal electrical resistivity. 17 

Observe that this study may provide a means of probing the phonon-mechanism in exotic 
superconductors, such as, heavy fermion superconductors, organic superconductors, fullerene 
superconductors, and high T c cuprates. For superconductors caused by the electron-phonon 
interaction we expect the following behavior. As the electrons are weakly localized by 
impurities or radiation damage, the electron-phonon interaction is weakened. As a result, 
both T c and TCR are decreasing at the same rate. When A is approaching zero, both T c and 
TCR drops to zero almost simultaneously. When this happens we may say that the electron- 
phonon interaction is the origin of the pairing in the superconductors. This behavior was 
already confirmed in A15 superconductors 14-17 and Ternary superconductors. 71 More details 
will be published elsewhere. 

VI. CONCLUSION 

It is shown that weak localization decreases both the conductivity and the electron- 
phonon interaction at the same rate and thereby leads to the Mooij rule. As the residual 
resistivity is increasing due to weak localization, so the thermal electrical resistivity is de- 
creasing, producing the decrease of TCR. When the electron-phonon interaction is near 
zero, only the residual resistivity is left and therefore the negative TCR obtains. This 
study may provide a means of probing the phonon-mechanism in exotic superconductors, 
such as, heavy fermion superconductors, organic and fullerene superconductors, and high T c 
superconductors. 
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Table I. Comparison of Xt r and the McMillan's electron-phonon coupling constant A. 
Data are from Economou, Ref. 38 and Grimvall, Ref. 39. 



Metal 




A 


Metal 




A 


Li 


.40 


.41±.15 


Na 


.16 


.16±.04 


K 


.14 


.13±.03 


Rb 


.19 


.16±.04 


Cs 


.26 


.16±.06 


Mg 


.32 


.35±.04 


Zn 


.67 


.42±.05 


Cd 


.51 


.40±.05 


Al 


.41 


.43±.05 


Pb 


1.79 


1.55 


In 


.85 


.805 


Hg 


2.3 


1.6 


Cu 


.13 


.14±.03 


Ag 


.13 


.10±.04 


An 


.08 


.14±.05 


Nb 


1.11 


.9±.2 
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FIGURES 



FIG. 1. The temperature coefficient of resistance a versus resistivity for bulk alloys (+), thin 
films (•), and amorphous (X) alloys. Data are from Mooij, Ref. 19. 

FIG. 2. Resistivity versus temperature for Ti and TiAl alloys containing 0, 3, 6, 11, and 33% 
Al. Data are from Mooij, Ref. 19. 

FIG. 3. McMillan's coupling constant A versus dp/dT oc X tr for Ag-Ga, Au-Al, and Au-Ga 
alloys. Data are from Rapp, Ref. 52 and Grimvall, Ref. 32. 

FIG. 4. (a) Phonon-limited resistivity p ph versus T for k F t = 15, 5, 3.4, 2.8, 2.4, and 2.2. (b) 
residual resistivity po versus T for the same six values of kpi. 

FIG. 5. Calculated resistivity versus temperature for fcp£ = 15, 5, 3.4, 2.8, 2.5, and 2.3. The 
solid lines are p(T) from an accurate formula, Eq. (31). The dashed lines represent the resistivity 
obtained from the approximate expression, Eq. (24). 
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